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Abstract. Discrete dynamical systems denned by the iteration of a 
polynomial map of the unit simplex to itself appear in the context of 
multilocus systems evolving under mutation, recombination and weak 
selection. Although remarkable progress has been made in finding par- 
ticular solutions to these systems, our knowledge of the general proper- 
ties of the space of all possible dynamical systems of this kind is very 
small. We prove that the space of bounded-degree polynomial maps of 
the unit simplex to itself is a compact and convex subset of a Euclidean 
space. We provide an explicit characterization of such a space and of its 
boundary. A special class of maps in the boundary, the folding maps, 
which generalize the logistic map for any dimension and degree are de- 
fined and constructed. Finally, we use numerical methods to study the 
ergodic and mixing properties of maps in the neighborhood of several of 
these folding maps. 



1. Introduction 

The dynamics of completely linked multilocus systems evolving neutrally 
under recurrent mutations are elegantly described by means of finite Markov 
chains. If the number of polymorphic genotypes in the population is finite 
and equal to n + 1, we can enumerate each genotype by an integer number i 
and denote by Xi the proportion of individuals that share this genotype. The 
genotype composition of the population {xi}™^ 1 is then described by a point 
in the n-dimensional unit simplex, or equivalently by a point in the space of 
probability distributions on the set of n + 1 genotypes. Assuming that the 
population mates at random with no genetic drift, the corresponding time 
evolution is simply given by 

n+l 

Xi (t + 1) = Y,Mi^ jXj (t), 
i=i 

where -Mj<_,- denotes the proportion of individuals of type j that mutate to 
type i in one generation. For this class of dynamical systems, the general 
structure of the space of evolutionary operators can be determined easily. 
In particular, as Y17=i ^i^j = 1 f° r eacn 3, an d every matrix element 
Mi<_j is non- negative, the space of Markov matrices is the direct product 
of n + 1 simplices of dimension n. We denote such a n(n + l)-polytope by 
PMaps(A n , 1), and by A n we denote the n-dimensional unit simplex. As the 
largest eigenvalue of M is 1 and because of the Perron-Frobenius theorem, 
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the dynamics of any Markov chain in the interior of PMaps(A n , 1) always 
converges to a unique fixed point in the n-simplex. This is because while 
one eigenvalue is 1, the remaining eigenvalues are in the interior of the unit 
disk in the complex plane. This compares with the Markov chains that are 
located in the boundary of the polytope PMaps(A n , 1), where the particu- 
lar arrangements of zeroes in the stochastic matrix span different faces. In 
particular, stochastic matrices located in faces of decreasing dimensionality 
correspond to matrices with an increasing number of zeroes. The spectrum 
of these matrices consists of one eigenvalue equal to 1, a set of eigenvalues 
located in the interior of the unit disk in the complex plane, and another set 
of eigenvalues located in the boundary of the unit disk. Indeed, when M^j 
has the maximum number of zero entries and its determinant is not zero, i.e. 
the matrix is a permutation matrix, all the eigenvalues are located in the 
boundary of the unit disk. The spectrum of any of these stochastic matri- 
ces allows us to determine some dynamical properties of the corresponding 
Markov chains; e.g. the existence of cycles is associated with more than one 
eigenvalue in the unit circle of the complex plane (see for instance |Sal97j ). 

In this note we investigate some analogous properties of the space of 
bounded-degree stochastic polynomial maps from the simplex into itself. 
First, we demonstrate that these spaces are compact and convex subsets 
of Euclidean spaces. This property is trivially true in the case of degree 
one stochastic maps, as we showed above. Second, we characterize a class of 
maps in the boundary of these spaces that generalize the notion of stochastic 
matrix with maximum number of zeroes and nonzero determinant. Geomet- 
rically, these polynomial maps are defined by maximizing the complexity of 
the algebraic subset in the preimage simplex that is mapped to the whole 
boundary of the simplex. Indeed, this is satisfied by a family of maps that 
we introduce here, the minimal polynomial folding maps, which, as we show 
below, are related to the theory of simplicial tessellations of the Euclidean 
space and to the affine Weyl groups [HW88, Vcs9lj. Finally, we show how 
these maps can be constructed by solving certain systems of polynomial 
equations, we compute several examples in one and two dimensions, and we 
study the dynamical properties of these maps and of their deformations in 
the space of polynomial maps. 

1.1. Basic notation and definitions. We denote by A" C l n+1 the stan- 
dard n-simplex: 
(1.1) 

A n = jz = (xi,x 2 ,...,£„+i) G M n+1 : J^Xi = 1, xj > 0, 1 < j < n + lj 

For simplicity, the vertices, edges, faces and other boundary components of 
the n-simplex are denoted as its 0-faces, 1-faces, 2-faces, . . . , and (n — 1)- 
faces. The interior of the n-simplex, A n \dA n , is also referred to as the 
n-face of the simplex. All the i-faces except for the 0-faces, which are closed 
sets in M n+1 , are open sets. Hence, any simplex can be expressed as the 
disjoint union of its faces. This implies that if S is the set of all the faces in 
A n , there exists an onto map F c : A" — > S, the face map, which assigns to 
each point in the simplex its corresponding face, i.e., x € F c {x). 



ON THE SPACE OF EVOLUTIONARY OPERATORS AND FOLDING DYNAMICS 3 

Any polynomial map from the n-simplex to itself is defined by means of 
n + 1 polynomials {Pi{xi,X2, • • • , Xn+i)}™^ satisfying 

n+l 

p i( x ) = 1 and p j( x ) > 0, 1 < i < n + 1, 

i=l 

for all x G A n . We define the degree of a given map as the degree of the 
polynomial Pi{x) of the largest degree. This definition is ambiguous because 
any pair of polynomials, Pi{x) and Pi(x)(x\ +X2 + • • • + x n+ \) N restricted to 
A n , are identical for any non- negative integer N . To avoid this ambiguity, 
we assume throughout that the degree of a polynomial P(x) is the smallest 
degree possible after factoring out (x\ + X2 + • • • + x Tl+ i) terms. We denote 
by PMaps(A n , k) the space of polynomial maps of degree less than or equal 
to k from A" to itself. By we denote the vector space of polynomials in 
n variables of degree at most k. 

A d- folding map from the simplex to itself is a d-to-1 and onto continuous 
mapping in which the number of preimages is d everywhere in the interior 
of A n . A polynomial fold of A n is a particular instance of a folding map 
that is also a polynomial map. 

2. The geometry of the space of stochastic polynomial maps 

For a given positive integer d, the corresponding d-folding polynomial 
maps located in the space of polynomial maps PMaps(A n , k) with the low- 
est degree k possible, are a very special set of points in PMaps(A n , k). 
This will become clearer after we describe in more detail the geometry of 
PMaps(A n ,/c). 

Theorem 1. The space of polynomial maps of degree k from the n-simplex 
to itself, PMaps(A n , k), is a compact and convex subset of the vector space 
of polynomial maps of degree k from R n to itself. 

Proof. The proof consists of an explicit construction of PMaps(A n , k), 
such that the properties stated in the theorem follow from such construction. 
For simplicity, we work with the projection of the simplex A n C M. n+1 down 
to R n , where A n is defined by the set of inequalities 

x% > 0, 1 < i < n, 

n 
i=l 

Any polynomial map / in PMaps(A n , k) can be defined by means of n poly- 
nomials {Pi(x)}2 = i of degree less than or equal to k, that obey the following 
conditions: 

(2.1) Pi{x) > 0, Vi G [l,n], Vx g A n , 

n 

(2.2) 1-J2 p i( x ) >°> VxG A n , 

i=i 

and 



(2.3) 



x' = f{x), Xi = Pi(x), 
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which states how preimages x are mapped to images x' and therefore, it 
completes the definition of /. The condition in Eq. (|2.ip can be simpli- 
fied thanks to an important theorem of Polya, which provides a systematic 
process for deciding whether a given polynomial is strictly positive on the 
simplex [HLP52] . Using Polya's result, one can show how the polynomials 
that are non- negative on A n span a convex cone in IF^, the space of poly- 
nomials in n variables of degree at most k. To show this, we first need to 
construct homogeneous representatives of the polynomials in 11^. In partic- 
ular, if P(x) is a polynomial in IljJ, we construct its homogenous polynomial 
representative Pr{x) in R n+1 as follows. We write P(x) as a sum of mono- 
mials 

\a\<k 

with \a\ = Y17=i a *' an< ^ multiply the terms with \a\ < k by {x\ + x<i + • • • + 
x n +i) fc_ ' a L In this way, we obtain a homogenous polynomial of degree k 
which is identical to P{x) when restricted to A" C M n+1 : Ph(%) = 

J2 C *i»2-an (*|«|,fc + (1 " $\a\,k)(xi +X2--- + X„+i)*- |a| ) X^ 1 X% 2 •••<", 
\a\<k 

with 5\ a \ t k = 1 when \a\ = k and S^ j, = otherwise. Using this homogenous 
polynomial Pjj(x), Polya's theorem simply says that if P(x) is positive on 
A n , i.e. P(x) > 0, Vx € A", then there exists a non-negative integer N 
such that all the coefficients of (xi + X2 + ■ ■ ■ + x n+ i) N Pjj{x) are positive. 
If P(x) is positive on A n then XP(x), with A € M + , is also positive. This 
implies that the polynomials of degree k with P{x) > 0, Vx € A™, span an 
open convex cone /C+(A n , k) C II£. It also follows that dim /C+(A n ,/c) = 
dim II?. Furthermore, one can use this construction to describe the set of 
non-negative polynomials on the simplex. More precisely, if P(x) is strictly 
positive on the n-simplex and Q(x) is non-negative for all x £ A n , then any 
deformation of Q(x) by eP(x) makes Q(x) + eP(x) an element of /C+(A n , k) 
for any positive real value of e. Therefore, the closure of lC+(A n ,k) in IF^, 
which we denote by JC(A n ,k) = cl (/C + (A n , k)), is the set of non-negative 
polynomials on A n . In other words, 

{P(x) € II£: P{x) > 0, Vx € A™} = K(A n ,k). 

Thus, the set of polynomials {Pj(x)}™ =1 that satisfy Eq. (|2.1|) are elements 
of the cartesian product of n closed convex cones lC(A n ,k). This product 
of cones is itself a closed convex cone in ©™ =i U]!, the direct sum of vector 
spaces of polynomials on the simplex. We denote such a convex cone by 
xf =1 /C(A",&). 

The space PMaps(A n , k) C ®™ =i ITjJ and its compactness property arise 
after we impose Eq. (|2.2p on the convex cone x™ =1 /C(A n , k). By the 
Heine-Borel theorem it is sufficient to show that PMaps(A", k) is a bounded 
subset of the cone x" =1 /C(A n , k) C ffi" =i nj?, in order to demonstrate that 
PMaps(A n , k) is compact. First, let 

span A {AP;(x)}™ =1 

for A € M + positive and real, denote the ray in x^ =1 lC(A n ,k) generated 
by {.Pj(x)}" =1 . The maximum of the function Y^i=x Pi( x ) on the simplex is 



ON THE SPACE OF EVOLUTIONARY OPERATORS AND FOLDING DYNAMICS 5 



strictly positive for any non-zero set of polynomials {Pi(x)}^ =1 £ x™ =1 /C(A n , k). 
Let us denote such a maximum value by S max . Now, Eq. (|2.2p restricted 
on this ray is equivalent to requiring the non- negativity of 1 — A^" =1 Pi(x) 
for all x in the simplex. This is only satisfied by the segment {\Pi(x)}™ =1 
parametrized by A € [0,1/ S max ]. Similarly, for every ray in x™ =1 /C(A n ,fc) 
only a bounded and finite subinterval of the ray satisfies Eq. (|2.2p . It fol- 
lows that PMaps(A n ,fc) is a bounded and closed subset of the convex cone 
xf =1 }C(A n ,k). 

The convexity of PMaps(A n , k) follows from the convexity of the simplex 
A n . In particular, if the image of x defined by {Pi(x)}f =1 G PMaps(A n , k) is 
in A n for each x € A n , and the same is true for a second map {Qi(x)}f =1 € 
PMaps(A n ,/c), then {tPi(x) + (1 - t)Qi(x)}^ =1 is in A n and yields a well- 
defined map in PMaps(A n , k) for any t G [0, !].□ 




FIGURE 1. Convex set of quadratic maps of the type f(x) = 
x{A+Bx) on the 1-simplex. The horizontal axis parametrizes 
the A-coefficient and the vertical axis parametrizes the B- 
coefficient. For pairs {^4, B} outside the convex set, the im- 
age of the quadratic map is not in the simplex. The red dot 
denotes the identity map, while the black dot denotes the 

! i i r ij- £ I \ 1— cos(2 arccos(l— 2x)) 

polynomial two-foldmg map j2\x) = 5 — = 

Ax{l - x). 

Our construction of PMaps(A n , k) relies on Polya's criterion to decide 
whether a given polynomial vanishes somewhere in the simplex A n or not. 
In particular, we have showed how the interior of PMaps(A n , k) consists of 
maps for which every defining polynomial, {Pi{x)}™ =l plus the polynomial 
1 — ^ ■ Pj(x), are strictly positive on A n , and therefore, their corresponding 
images are in the interior of A n . This compares with the boundary of 
PMaps(A n , k), consisting of maps for which some of the defining polynomials 
do vanish on some algebraic subsets of the standard simplex. For these 
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maps the preimages of the boundary of A n are not the empty set. The case 
of finite Markov chains is particularly illuminating in this regard. When 
k = 1, the defining polynomials are non- negative linear forms on the simplex, 
while the boundary of PMaps(A n , 1) consists of maps that have a non-empty 
preimage of the boundary dA n . In an extreme limiting case, there exist 
maps with non-singular Jacobian matrices that maximize the size of the 
corresponding preimages of dA n . These maps consist of invertible stochastic 
matrices that have the largest number possible of zero entries, i.e. the 
permutation matrices. Geometrically, these matrices are special vertices in 
the boundary of the n(n + l)-dimensional polytope PMaps(A n , 1), among 
other things because they are the only bijective maps therein (see Lemma 1 
in the Appendix). 

One can define nonlinear analogs of these extremal maps in PMaps(A n , k) 
for k > 1. They consist of polynomial maps /: A n — > A n with a non- 
singular Jacobian almost everywhere in A n , that maximize the complexity 
of the algebraic subset in A" that is the preimage of the boundary under 
/, i.e. they maximize the complexity of f^ 1 (<9A n ) C A n . These maps are 
easy to construct in the case of the 1-simplex. For instance, when n = 1 
and k = 2 the largest algebraic subset of the interval that can be mapped 
to its boundary by a quadratic map (with a non-singular Jacobian almost 
everywhere in [0,1]) consists of three points (see Fig. 1). In general, the 
largest algebraic subset of A 1 that can be mapped to the boundary of A 1 
by a map in PMaps(A 1 ,/c) consists of k + 1 points. These maps can be 
constructed explicitly [HW881 IVes91j . They coincide with the Chebyshev 
maps and their permutations: 



Here, T)% denotes the k Chebyshev polynomial of the first kind, which can 
be evaluated as T k (x) = cos (fearccosx). One can show from the properties 
of the Chebyshev maps that every map fk is a k- folding map of the 1-simplex. 

As far as we know, there do not exist previous attempts in the literature 
to construct the higher dimensional analogs of these extremal maps on the 
unit simplex. Although these maps are indeed harder to construct, we show 
below how some particular families of extremal maps admit a simple geo- 
metrical description. In particular, the concept of a polynomial map / of 
fixed degree that maximizes the complexity of / _1 (dA n ) finds an instance 
in any polynomial (i-folding map of minimal degree. As every polynomial 
d-folding map / is continuous, onto and d-to-1 everywhere in the interior 
of the image simplex, one can decompose the preimage simplex A n as the 
union of d sets, A™ = U^ =1 5 a , such that the restriction of the map on each 
subset, f\s a - 5 a — > A n , is one-to-one and onto. In this case, the preim- 
age of the boundary is simply the union of the boundaries of each subset 
$a, f~ l (<9A n ) = U^ =1 d5 a . Therefore, those d-folding maps / located in 
spaces PMaps(A n , k) that have the smallest k possible maximize the com- 
plexity of f~ l (<9A n ) C A n among all maps in PMaps(A n , k). In the next 
section we define and construct these minimal polynomial folding maps. 
First, we propose a restricted notion of homotopy that allows us to group 
together equivalent folding maps. For instance, we consider 4x(l — x) and 



(2.4) 



f k : [0,1] -> [0,1], f k (x) 



l-r fc (l-2a) 
2 
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2x0(1/2 — x) + (2 — 2x)0(x — 1/2) to be topologically equivalent 2- folding 
maps of the 1-simplex (here, denotes the Heaviside function); however, 
we do not consider the polynomial 2-fold 4a; (1 — x) and the piecewise linear 
3-fold 3x0(1/3 - x) + (2 - 3x)0(x - 1/3)0(2/3 - x) + (3x - 2)0(x - 2/3) to 
be equivalent. Then, we use this notion of homotopy to define the minimal 
polynomial folding maps. Finally, we show how to construct these maps as 
the zero loci of certain polynomial systems equations and compute several 
examples. 

3. Polynomial Folding Maps 

In order to be able to compare polynomial folds with non-polynomial 
folds, e.g. piecewise linear folding maps, we need an appropriate criterion 
to establish a topological equivalence relation between them. To that end, 
we introduce a restricted notion of homotopy equivalence between folding 
maps as follows. We say that two maps, / : A n — > A n and g : A n — > A n , are 




Figure 2. A piecewise linear 2-fold of the 1-simplex. 
Here, the 1-simplex spans the interval [0,1], the horizon- 
tal axis denotes the preimage of the map, and the verti- 
cal axis denotes its image. The standard notion of ho- 
motopy equivalence between maps is not valid to compare 
folding maps. The 2-folding map, given by the function 
f(x) = 1Ox0(^-x) + (1-^(x-^))0(x-^), is homotopic 
to the polynomial fold g{x) = 4x(l — x) (see left plot). The 
same 2-folding map is also homotopic to the identity map 
I{x) = x, (see right plot). 

face- constrained homotopic (or fc-homotopic) if and only if there exist two 
continuous functions, H : A n x [0, 1] — > A n and K: A n x [0, 1] — > A n , from 
the product of the n-simplex with the unit interval [0, 1] to the n-simplex 
such that: 

i) K(x, t) : A n — > A n is one-to-one and onto for all t £ [0, 1]. 
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ii) {K(x,0), H(x,0)} = {x, f(x)} and {K(x,l), H(x,l)} = {y, g(y)} 
for all x e A n . 

iii) The face- constrained condition requires that if F c (x) is the face 
where x belongs, x G F c (x) C A n , and F c (f(x)) denotes the face 
where the image f(x) is located, f{x) G F c (f(x)) C A n , then the 
deformation defined by K and H preserves such constraints, i.e. 
K(x,t) G F c (x) and H(x,t) G F c (f(x)) for all t G [0, 1]. 

Note that if the condition iii) is not satisfied, one recovers the standard 
notion of homotopy by choosing K to be the identity map. 




Figure 3. A face- constrained homotopy of the piecewise 
linear fold defined by f(x) = WxQ — x) + (1 — ^(x — 
jq))Q(x — jq). The preimage points x = and x = 1 are 
located on 0-faces, while both of their images are on the same 
0-face (y = 0). The images of all the interior points, except 
x = 1/10, are also interior points. The image of the interior 
point x = 1/10 is located on the 0-face y = 1. These con- 
straints are preserved by any fc-homotopy. Hence, while f(x) 
is fc-homotopic to g(x) = 4x(l — x), it is not fc-homotopic to 
the identity map. 

Definition 1. Minimal polynomial folding map. Given an fc-homotopy class 
of folding maps on the n-simplex, we say that a particular polynomial repre- 
sentative is minimal if and only if the corresponding polynomial map exists 
in PMaps(A n , k) and there do not exist fc-homotopically equivalent maps in 
PMaps(A n ,/), with I < k. 

Let us assume that a given minimal polynomial folding map / : A n — > A n 
is defined by the polynomials {Pj(x)}" =1 and P n+ \{x) = 1 — Y^j=i Pj( x )- 
Then, the preimage of each of the n + 1 facets of the simplex is an algebraic 
subset of the unit simplex defined either by 

{x G A n : Pi(x) = 0}, 
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for the facet = 0, or by 

n 

{x G A n : 1 ~Y, P j( x ) = p n+i(x) = 0}, 

J'=l 

for the facet Y%=iVi = 1- We denote by {z; C A n \dA n , 6./ C dA n } the 
co dimension-one preimage of the i th facet, decomposed in its interior and 
boundary parts. In this case, Zi is the subset of the preimage located in the 
interior of the simplex and bi is the complementary subset located in the 
boundary of the simplex. It is important to note that as 6j C dA n is the zero 
locus of a polynomial defined on A n , bi is either the empty set, or a finite set 
of facets. Using other words, bi cannot be composed of finite subsets within 
facets because such subsets are not algebraic; bi can only be composed of 
whole facets, because they are the only algebraic subsets of codimension 
one in the boundary of the simplex, e.g. Xi = 0, or 1 — X^j=i x j = 0- 
Indeed, there exist folding maps that do not admit fc-homotopic polynomial 
representatives; this is because whenever different subsets within a single 
facet are mapped to distinct facets, the folding map cannot be polynomial. 
This will become clearer after we show the following. 

Theorem 2. Let f £ PMaps(A n , k) be a minimal polynomial folding map 
of degree k. Then, each defining polynomial Pi(x), with 1 < i < n+ 1, can 
be factored as 

(3.1) Pi(x) = li(x)(qi(x)) 2 mi(x), 

where U{x), qi(x) and rrii(x) obey the following: 

• The zero locus of U{x) in A n is exactly bi C dA n . 

• The degree of li (x) equals the number of facets contained in bi . Fur- 
thermore, the product of all the U{x) 's is 

n+l n 
i=l j i=l 

• The zero locus of qi{x) in the interior of the n-simplex is exactly 
Zi C A n \9A n . 

• U(x) and mi(x) are strictly positive in the interior of the n-simplex. 
The zero set of m,i(x) is either empty or it is contained in the zero 
set of li{x). 

• The degree of Pi(x) is bounded by k, deg/j + 2deggj + deg rrii < k. 

Proof. The preimage of A n \<9A n under / consists of d disconnected open 
subsets in the interior of A n . The restriction of / to any of these subsets 
gives rise to a one-to-one mapping that has the whole interior of the simplex 
as codomain. The complementary set of the union of these d open subsets in 
A n is f~ 1 (dA n ). These facts follow from the properties of the folding maps: 
they are continuous and onto, for every point in the interior of the simplex 
there exist d different preimages, and they map the boundary <9A n to a 
subset in the boundary of the image simplex (see Lemma 2 in the Appendix). 
In addition to this, the fact that a minimal polynomial folding map is also 
polynomial imposes further constraints. In particular, the properties of the 
defining polynomials stated in Theorem 2 are derived from the algebraic 
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properties of the preimage set / _1 (<9A n ). Algebraically, f~ l (<9A n ) is the 
union of irreducible varieties of codimension one in A n . The set defined by 
the zero locus Pi(x) = denotes the algebraic variety in f" 1 (<9A n ) that 
is preimage of the i th facet. As such a preimage has to be the union of a 
finite number of facets and a codimension one variety in the interior of the 
simplex, Pi{x) can be factored as 

Pi(x) = li(x)n(x)mi(x). 

Here, the zero locus of U{x) is the union of facets in the boundary, the zero 
locus of ri(x) is a codimension one variety in the interior of the simplex, and 
rrii{x) is a non- negative polynomial of degree equal or less than k — deg li{x) — 
degrj(x), whose zero locus is either empty or it is contained in the zero locus 
of k(x). The existence of rrii(x) G /C(A n , k — deg^(x) — degr,(x)) follows 
from the facts that Pi{x) has a degree equal or less than k, that Pi(x) > 
for all x € A n and that the zero set of Pi(x) is already specified by r,(x) and 
li(x). Because of these constraints, the existence of an additional polynomial 
factor of degree less than or equal to k — degij(x) — degrj(x) is possible. 

As the i^-facet of the n-simplex is defined by the zero locus of the degree- 
one monomial Xi (when i = n + 1, x n+ \ = 1 — Yll=i the l° wes t degree 
polynomial whose zero locus consists of a predefined set of facets 6j C A" is 

k(x) = Xj. 

If bi is empty then U{x) = 1. It follows that U(x) > for all x in the interior 
of the simplex. Furthermore, as f(dA n ) is contained in <9A n , the images of 
every facet in A" are in the union of all the image facets, and therefore 

n^) 

i=l 

has to be a polynomial divisible by (1 — X^j=i x j) Yli=i x i- As the preimage 
of each facet is the union of irreducible varieties of codimension one, each 
facet has to be mapped to a whole facet. This implies that the monomial 
Xj whose zero locus defines the facet cannot appear as a factor in more 
than one polynomial U{x). Otherwise, the image of such j^-facet would 
be a face of codimension larger than one. It follows that k( x ) is 

not only divisible by (1 — Y^j=i x j) YYl=i x i but it i s a constant multiple of 
(1 — Sj=i x j) lir=i x i- By convention, such a constant factor is absorbed in 
the mi(x) terms, and 

ra+l n n 

i=l j=l 1=1 

Finally, it remains to show that the factor ri(x) is the square of a poly- 
nomial: ri(x) = (qi(x)) . As the zero set Zi = {x E A": ri{x) = 0} is a 
co dimension-one variety in the interior of the simplex, and rj(cc) is strictly 
positive for every x € A n \zi, it follows from the Taylor expansion of r^(x) 
around any x € Z{ that the gradient of ri{x) restricted to Z{ has to vanish. 
Therefore, Zi has to be the zero locus of a polynomial q%(x) which might have 
negative values on a subset of the simplex and satisfies ri{x) = (qi(x)) 2 .d 
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A priori, the factorization in Eq. (|3.ip can be used to construct any 
minimal polynomial folding map. First, one fixes an fc-homotopy class of 
folding maps that admits polynomial representatives. Then, the assignment 
of facets to image facets defined by the fc-homotopy class is used to deter- 
mine the polynomials {^0c)}2=a ■ Unfortunately, specifying the remaining 
factors in Eq. (|3.ip is a harder problem. The approach that we follow con- 
sists of parametrizing the set of all the polynomials, {rrii(x, (f) l )}1=i and 
{qi(x, 0)}™=i, that are compatible with the given class of folding maps. 
Then, we determine the particular parameters {6, (p 1 }™^ that define the 
folding map by solving a system of polynomial equations. More precisely, 
we use the defining equation 
(3.2) 

n 

l n +i{x) (q n+ i(x,6)) 2 m n+1 (x,(j) n+1 ) = 1 - } )k(x) (qi{x , 0)) 2 rrii(x , <j?) , 

to determine the parameters {8, (i> % Yi=\ that satisfy Eq. (|3.2|) for all x € A". 
In practice, this is a tedious approach because the set of compatible polyno- 
mials usually has several independent components and each component is a 
complicated semialgebraic set in a high-dimensional vector space. Neverthe- 
less, it is possible to identify each component and construct the correspond- 
ing spaces of polynomials. In particular, every component is associated with 
a different partition {(aj,6j)}j of the degrees of (qi(x)) 2 mi(x); i.e. partitions 
of k — deg li(x) = 2di + bi, with ai = deg qi(x), bi = deg rrii(x) and k the de- 
gree of the polynomial map. Given a particular partition {(oi,6i), (02,62), 
. . . (a n +i, b n+ \)} of the degree of every polynomial factor, it is possible to 
construct the space of compatible polynomial factors with these prescribed 
degrees and find the solutions of Eq. (|3.2I) . 

This construction becomes more transparent in the next subsection, where 
we show several examples in one and two dimensions. However, before we 
construct particular examples, it is useful to state some general properties 
of these families of polynomials. For instance, the family of polynomials 
{qi(x, 0)}™^ can be determined by requiring the corresponding zero loci 
{z{ C A* 1 }™^ 1 to be compatible with the given class of folding maps. This 
allows for the possibility of the existence of maps in the fc-homotopy class for 
which the preimage of dA n restricted to the interior of the simplex is exactly 
defined by the set of hypersurfaces {zj}™^ 1 . The parameter space {<f> 1 } 
admits a simpler description. The set of compatible polynomials m,j(x,0*) 
consists of those polynomials of degree bi 

m i(x,4>)= ( Pa 1 a2---a n x l lx 2 2 ' ' ' X n" 

\a\<bi 

whose zero loci are either non-existent or are contained in the zero locus of 
k(x). 

3.1. Examples of minimal polynomial folding maps. In the following 
examples we fix the fc-homotopy class where the minimal polynomial folding 
map sits by defining a non-polynomial folding map in the same class of 
maps. In particular, a well studied family of folding maps are the piecewise 
linear folding maps derived from periodic tessellations of the n-dimensional 
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Euclidean space |HW88j . In this construction one defines an embedding of 
the simplex into a Euclidean space of the same dimension, such that the 
action of the group of reflections through the facets of the simplex generates 
a tessellation of the space (e.g. see Fig. 4). One then uses this tessellation 
to determine the associated folding maps. More precisely, we can construct 
these piecewise linear folding maps by means of (see |HW88j for more details 
on this construction): 

i) An affine embedding of A n in M. n such that the embedded simplex 
tessellates M n by means of reflections through the hyperplanes asso- 
ciated with the (n — l)-faces of the simplex. We denote this affine 
map by A: A n — > R n , and we assume that the image of one of the 
vertices of A n is the zero vector in ~R n . 

ii) An integer dilation of the image simplex in M. n . In particular, we 
require that for any integer dilation, the dilated simplex is tessellated 
by a finite number of smaller replicas of itself (e.g. see Fig. 4.). 

iii) A map that assigns to each point in every tile of the tessellation of 
~R n the equivalent point in the tile that is the image of A n under A. 
This map is unique; it is constructed by means of reflections through 
hyperplanes; and we denote it by h: W 1 — > A(A n ). 

Now, we define a folding map by first dilating the affine simplex A(A n ) m 
units, i.e. we multiply every vector in ^4(A n ) C W 1 by the positive integer 
m. We denote such dilation by m o A(A n ), and by d we denote the number 
of smaller replicas of the simplex contained in m o A{A n ). Therefore, the 
piecewise linear d-folding map is defined as the following composition of 
mappings 

A' 1 o homo A: A n -> A n . 

In what follows we define several piecewise linear folding maps in one and 
two dimensions, and then apply the polynomial factorization of Theorem 2 
to determine their associated minimal polynomial folding maps. 

3.1.1. Folding maps of the one-simplex. The piecewise linear (i-folding maps 
in dimension one can be derived from partitions of the unit interval in d 
subsegments. In particular, if we assume that x = is mapped to then a 
partition determines uniquely a piecewise linear folding map. The simplest 
non-trivial case, the two-folding map, can be written as f-2 a {x) = (x/a)Q(a— 
x) + ((1 — x)/(l — a))Q(x — a) with a G (0, 1); although for simplicity we 
consider only a = 1/2. It follows from the definition of /2, a that there does 
not exist a strictly linear map that is fc-homotopic to this two- folding map. 
However, if there exists a polynomial map of degree two that is fc-homotopic 
to the two-fold, then by Theorem 2 its corresponding defining polynomials 
will factor as 

(3.3) P 1 {x)=Ax{l-x), P 2 {x) = l-P 1 {x) = {B + Cx) 2 , 

with A > 0, B and C being unknowns to be determined. The fc-homotopy 
equivalence to f2,a{ x ) requires that the preimages of y = are x = and 
x = 1, and that there exists an interior point Xj = —B/C G (0,1) that 
is mapped to y = 1. We can determine the coefficients A, B and C by 
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expanding the denning equation P2(x) = 1 — Pi{x) in powers of x 

(B + Cxf = 1 - Ax(l -x), B 2 - 1 + {2BC + A)x + (C 2 - A)x 2 = 0. 

If this is true for every x 6 [0, 1], then B = 1, C = -A/2, A 2 /A = A, and 
A = 4. Therefore, the corresponding minimal polynomial folding map is 
defined by P\(x) = 4x(l — x) and P^,{x) = (1 — 2x) 2 , i.e. the logistic map 
|May76| . 



k 


Pi(x) 


P 2 (x) = l-Pi(x) 


1 


X 


1 — X 


2 


4x(l — x) 


(1 - 2x)^ 


3 


x(3-4x) 2 


(1 -x)(l -Ax) 1 


4 


16x(l -x)(l - 2x) 2 


(8x 2 -8x + l) 2 


5 


x(5 - 20x + 16x 2 ) 2 


(1 -x) (I6x 2 - 12x + l) 2 


6 


4x(l-x)((l-4x)(3-4x)) 2 


((l-2x)(l-16x + 16x 2 )) 2 



Table 1. Polynomial factorization of the first six Cheby- 
shev maps from the 1-simplex to itself. 



One can repeat this construction for every fc-homotopy class of d-folding 
maps on the one-simplex. As an example, we computed the first six mini- 
mal polynomial folding maps (see Table 1). It is easy to notice that these 
polynomials are the normalized Chebyshev polynomials of the first kind 
Pi,d(x) = 1 ~ co& ( dari ^ cos ( 1 ~ 2x )) , This is indeed expected from the definition 
of the Chebyshev polynomials, which employs the piecewise linear folding 
function [AR64] . 

3.1.2. Folding maps of the two-simplex. From the theory of affine Weyl 
groups, we know that every tessellation of the two-dimensional plane that 
is derived from the reflections of an isosceles triangle is isomorphic to one 
of two possibilities (see Fig. 4.). We denote each possibility by the name of 
its corresponding affine Weyl group (A2 and B2). Of these, only the fold- 
ing maps derived from the B2 tessellation admit polynomial representatives. 
This is because the folding maps derived from the A2 tessellation map every 
edge in the boundary of the triangle to more than one different edge, and 
therefore the pre- image of an edge cannot be an algebraic set. 

Given the fc-homotopy class associated with the B2 two-fold, using The- 
orem 2, and assuming that the minimal polynomial folding map has degree 
k = 2, the defining polynomials have to be of the form 

/ Pi(x,y) \ / A(y-Bx) 2 \ 

(3.4) P 2 {x,y) = (l-x-y)(C + Dx + Ey) . 

\l-P 1 (x,y)-P 2 (x,y) J \ Fxy J 

Here, A,B,...,F, are unknown parameters that satisfy A > 0, B > 0, 
F > and C + Dx + Ey > for all (x, y) € A 2 . Now, we impose Eq. ([33]) 
on this family of polynomials, which gives rise to the equation 

Fxy = 1 - A(y - Bx) 2 - (1 - x - y)(C + Dx + Ey). 
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Figure 4. Two tessellations of the real plane by 2- 
simplices. The tessellation on the left is associated with the 
Weyl group A2 and the equilateral triangle. The right picture 
represents a tessellation by isosceles right triangles, which is 
associated with the Weyl group B2. Both tessellations give 
rise to corresponding piecewise linear folding maps on the 
2-simplex. The dashed green triangles associated with the 
g = 1 dilation, give rise to the 1-fold; the dashed blue tri- 
angles correspond to g = 2 dilations associated with 4- folds; 
the red g = 3 dilations give rise to 9-folds. Of these folds only 
the ones associated with the isosceles right triangles admit 
fc-homotopic polynomial maps. 



Expanding this equation in the standard basis of polynomials of degree two 
with two variables, we get = 

{l-C) + {C-D)x+{C-E)y+{E+D+2AB-F)xy+{D-AB 2 )x 2 +{E-A)y 2 . 

Therefore, finding the parameters (A, B, . . . ,F) that satisfy this equation 
for all (x,y) € A 2 amounts to solving the following system of polynomial 
equations: 

1 - C = 0, 
C-D = 0, 

C-E = 0, 

1 ' E + D + 2AB-F = 0, 

D - AB 2 = 0, 
E-A = 0. 

From the first, second and third equations it follows that C = D = E = 1, 
which yields the strictly positive linear form 1 + x + y € /C_|_(A n , 1). From 
the last equation we find that A = 1, and therefore F and B have to obey 
the fourth and fifth equations with A = C = D = E = 1. This implies that 
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2 = F — 2B and 1 — B 2 = 0. As B has to be positive in order to yield a 
folding map compatible with the fc-homotopy class, B has to be 1 and F 
has to be 4. This shows that there is a unique polynomial two-fold of degree 
two in this fc-homotopy class. We can compactly write this map as 

/ P 1 (x,y) \ ( (y-x) 2 \ 

(3.6) P 2 (x,y) = (l-x-y)(l + x + y) . 

\l-P 1 (x,y)-P 2 (x,y) J \ Axy ) 



f d ={Pi(x,y), P 2 {x,y), P 3 {x,y) = l-P 1 (x,y)-P 2 (x,y)} 

fi = {x,y,l-x- y} 

h = {( x ~ y) 2 , (1 ~ x - y)(l + x + y), Axy} 
h = {(1 " - 2y 2 ) 2 , 8xy(l - 2xy), 

4(1 - x - y)(l + x + y)(x - y) 2 } 

/ 8 = {(l - 4x 2 + 4x 4 - 8xy - Ay 2 + 24x 2 y 2 + 4y 2 )' z , 
8(1 - x - y)(l + x + y)(l - 2x 2 + 2x 4 + 4xy - 2y 2 - 4x 2 y 2 + 2y A )(x - y) 2 , 

32xt/(1 - 2xy)(l - 2x 2 - 2y 2 ) 2 } 

f g = { x (l - y )(2-9x + 2Ax A - 16x a + 9y - 24y z + 16y a )(3 - 4x^(1 - 4yf, 
y(l - x){2 - 9y + 24y 2 - 16y 3 + 9x - 24x 2 + 16x 3 )(3 - 4y) 2 (l - 4x) 2 , 
(1 - x - y) 2 (l - 8x + I6x 2 -8y- 16xy + I6y 2 ) 2 } 

Table 2. Polynomial factorization of five different folding 
maps from the 2-simplex to itself. 



One can repeat this construction for higher degree folding maps (see Table 
2 and Fig. 5). The minimal polynomial folding maps of degree four and eight 
are unique and they can be expressed as a composition of the two-folding 
map; i.e. fi = f 2 ° ji and fs = f 2 o f 2 o f 2 . In the case of the nine- 
folding map we only explored the space of maps corresponding to a single 
partition of the degrees: {(2,4), (2,4), (2,4)}. In this partition the interior 
pre-images of the boundary correspond to quadratic curves (see last row in 
Table 2 and Fig. 5). Although we found a unique minimal folding map, 
other different partitions of the degrees in the defining polynomials could 
give rise to additional minimal maps. The fact that there may exist more 
than one equivalent minimal folding map can be exemplified by two folding 
maps of degree eighteen: fg o f 2 and f 2 o /g, which are at the same time 
distinct and fc-homotopically equivalent maps. 

3.1.3. Dynamical properties of some polynomial folding maps. The set of 
Chebyshev maps on the one-simplex has very rich dynamical properties 
[AR64j . First, one can show using simple trigonometry that every Cheby- 
shev map of degree d preserves the measure 

, „s 1 f dx „ 
fi(B) = - / BC 0,1. 

In addition to this, the corresponding dynamical systems are strongly mixing 
with respect to \i\ they form a commutative semigroup under the composi- 
tion of maps (fd°fe = fe°fd = fdxe)', an d they are orthogonal as polynomial 
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Figure 5. Several folding maps on the 2-simplex. The lines 
within each triangle denote the different preimages of the 
boundary of the simplex under a given folding map. Each row 
consists of an fc-homotopy class of folding maps represented 
by a piecewise linear map (see left) and a minimal polynomial 
folding map (see right and Table 2). The orange dots denote 
the fixed points of the minimal polynomial folding maps. 
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functions under the inner product L 2 ([0, 1], fx) denned by the invariant mea- 
sure fx |AR64j . Some of these properties are also present in deformations 
of the Chebyshev maps that are not folding maps. For instance, it is well 
known that many ergodic properties of the logistic map /2 are preserved by 
deformations of the type = Xx(l — x), where A belongs to a subset of 
full measure in the interval ~ (3.57,4] (see |May76| ). 

Despite these results in dimension one, our knowledge regarding polyno- 
mial folding maps on higher dimensional simplices is non-existent. Because 
of this, we performed several numerical experiments with the polynomial 
folding maps of degree two and nine that we derived above (Table 2). In the 
case of the polynomial two- fold, we determined the fixed points of the first 
three iterations of the map (i.e., fa, and fs)- We found that the spec- 
trum of the Jacobian matrix associated with these fixed points lies outside 
the unit disk in the complex plane, which is suggestive of chaotic behav- 
ior. We iterated 10 4 times the two-folding map for different random initial 
conditions, and found that the dynamics never converged to an attracting 
periodic orbit. We repeated this numerical analysis for thousands of maps 
in the neighborhood of /2 in PMaps(A 2 , 2), and found that the closest maps 
to /2 exhibited either ergodic behavior or had an attracting periodic orbit 
of size larger than 10 4 iterations (see Fig. 6). In order to perform these 
experiments, we endowed PMaps(A 2 , 2) with an 1? metric, and computed 
deformations of the two-folding map by constructing explicitly a polytopic 
approximation of the interior of PMaps(A 2 , 2) (see Appendix B for details). 

In the case of the nine-folding polynomial map we found thirteen fixed 
points. The Jacobian matrices associated with eleven of these fixed points 
have their spectrum located in the exterior of the unit disk, while in the 
remaining fixed points, (1, 0) and (0, 1), the Jacobian matrix has eigenvalues 
in the interior of the unit disk. These vertices are attracting fixed points, 
and each one is located next to two repelling fixed points (within a Euclidean 
distance of less than 1/100). We computed the dynamics of 10, 000 random 
initial conditions close to the origin (0, 0) and computed their fixation times 
(see Fig. 7). As expected, every scenario converged either to the vertex (0, 1) 
or (1,0), and the distribution of fixation times closely follows a log- normal 
distribution with mean 2.557 and standard deviation 0.482. 



4. Discussion 

Discrete-generation multilocus models with non-trivial linkage maps are 
difficult to study because of the nonlinearity associated with the dynami- 
cal equations. For instance, simple neutral mutation-recombination models 
with no genetic drift require quadratic maps from the simplex to itself to 
describe the corresponding population dynamics. More generally, the addi- 
tion of weak selection effects combined with time-dependent intensities of 
selection, mutation and recombination can give rise to very complex dy- 
namical systems also described by the iteration of polynomial maps from 
the simplex to itself (e.g. see |Lyub92j IG?o07l IBaaOll INHB99j ). Although 



these maps introduce a plethora of technical problems that are not present 
in the case of finite Markov chains, some particular solutions can be found 
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Distance to the two-folding map 



Figure 6. Dynamical properties of deformations of the 
polynomial two-folding map on the 2-simplex (see Appendix 
B). The horizontal axis denotes the L 2 distance to the two- 
folding map. The vertical axis denotes the smallest absolute 
value of the eigenvalues of the Jacobian matrix associated 
with the fixed points of the map. The green dots denote 
maps that are either ergodic or their dynamics converge to a 
periodic orbit of size larger than 10 4 iterations. The red dots 
denote maps whose dynamics converge either to a finite set 
of absorbing fixed points or to a periodic orbit of size smaller 
than 10 4 iterations. 



in biologically interesting regimes, such as the loose-linkage limit |NHB99j 
and the tight-linkage low- mutation limit [BaaOl]. 

In this paper, rather than stressing particular models of biological inter- 
est, we have described some general properties of the space of stochastic 
polynomial maps on the simplex. Our motivation here is closer in spirit to 
that of S. N. Bernstein in his work on stationary quadratic mappings from 
the simplex to itself |Be22} [Lyub92| . There, Bernstein was interested in the 
classification of quadratic maps Q £ PMaps(A n , 2) that reach equilibrium 
after one single iteration, QoQ = Q, generalizing the Hardy- Weinberg prin- 
ciple. Here, we have studied the global geometry of the space of polynomial 
maps PMaps(A n , k). We have introduced a systematic method to construct 
the interior of such a space of maps and characterized an important class of 
maps in its boundary, the minimal polynomial folding maps. The dynam- 
ics of some of these folding maps and their deformations in PMaps(A n , k) 
can exhibit ergodic and mixing properties, similar to the logistic map in di- 
mension one | May 76 1 . Therefore, any evolutionary model that falls in these 
regions of the space of stochastic maps would display striking dynamical 
properties with important consequences for the generation and maintenance 
of genetic variation in the population. As far as we know, our discovery 
of the existence of the minimal polynomial folding maps is new. We have 
demonstrated the existence of some of these folding maps by means of an 
original construction which demands significant computational resources. 
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Probability Density 




Fixation Time (generations) 



Figure 7. Distribution of fixation times associated with 
the dynamics of a minimal polynomial nine-fold (see Table 
2). Here, time is measured as the number of iterations. The 
vertex points (1,0) and (0,1) of the triangle (see Fig. 5) 
are the only stable fixed points of the map. The distribution 
represents the fixation times of 10 4 uniformly distributed ran- 
dom points in the square < x < 1/100, < y < 1/100. The 
shape of the resulting distribution did not change appreciably 
when smaller squares were chosen. A fit of the log-normal 
distribution (/t = 2.557, a = 0.482) to the simulated data 
(green curve) is shown superposed. 

Given that the one-dimensional minimal polynomial folding maps can be 
derived easily using classical methods, it remains an open question whether 
the higher dimensional folding maps can be derived using simpler methods 
than the ones that we use here. 
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Appendix A. Some Elementary Lemmas 

Lemma 1. Every one-to-one and onto stochastic matrix in PMaps(A n , 1) 
is an n + 1 by n + 1 permutation matrix. 

Proof. Let M be a one-to-one and onto stochastic matrix, represented by 
an n+1 by n+1 invertible matrix. Let the standard simplex A n be defined as 
the convex hull of the standard basis elements inM n+1 ,i.e. e x = (1,0,0,- ••), 
e2 = (0, 1, 0, • • • ), ... e n+ i = (0, 0, 0, • • • , 1). Therefore, the image of A n 
under M in M n+1 is the convex hull of the vertices {Mei}^ 1 . As M is 
a stochastic matrix, the image of the standard simplex is a subset of the 
standard simplex. Furthermore, as M is onto, the image simplex is indeed 
the standard simplex, i.e. the convex hull of the standard basis elements. 
This implies that Mei = ej for every i G [l,n + 1] and j £ [l,n + lj. As 
M is invertible, the image of the standard basis elements under M spans 
all the standard basis elements, and therefore M has to be a permutation 
matrix. □ 



Lemma 2. Let f : A n —¥ A n be a folding map. The image of the boundary 
of the simplex under f is contained in the boundary, f(dA n ) C dA n . 

Proof. It follows from the definition of folding map that / is continuous, 
onto and d-to-1, such that the number of preimages is d everywhere in the 
interior of A n . Let us assume that there exists a point x € dA n whose image 
under / is in the interior of the simplex. Then, one can construct a small 
open ball a ls° i n the interior of the simplex that contains fix). As 

fix) is the image of a point in the boundary of A n we can divide Bfr x ) in 
three subsets. First, one subset is an open disk of dimension n — 1 that is 
the image under / of a subset in the boundary of A n that contains x. This 
n — 1 dimensional disk, that we denote by Df( x )i divides the ball in two 
open half-balls B^ and BJ^, such that Bf^ = B^UD/^UB^. By 
the continuity of /, one of the half-balls B^ x ^ has a preimage open subset 
in A n whose boundary contains x G A n C M. n . However, the other half-ball 
BJ/ X \ does not have a preimage in A" whose boundary contains x. It follows 
that the number of preimages of x is larger than the number of preimages of 
points in Bj,y This is a contradiction because / is a folding map, and by 
definition the number of preimages is the same everywhere in the interior of 
the simplex. □ 



Appendix B. Deformations of maps 

In order to know whether the dynamical properties of a given folding map 
/* are preserved in its neighborhood in PMaps(A n , k), we need a systematic 
method to generate deformations of /* and determine their properties. First, 
we endow PMaps(A n , k) with the I? distance associated with the Euclidean 
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volume form on the simplex, i.e. 



9\\L* 



\ 



f^(P t (x) - Qi{x)f 



\dxi 



where /*, g € PMaps(A", k), and {Pi(x)}i and {Qi(x)}i are the sets of 
defining polynomials for /* and g respectively. Equipped with this metric 
we can generate deformations g S -B/„,e of the folding map within a distance 

B u = {g G PMaps(A n , k) : ||/* - g\\ L * < e}, 

by following the steps outlined in Fig. 8. In particular, we first generate 
random interior points r € PMaps(A n , k) and then take convex combinations 
of those points with the folding map to obtain the desired deformations 
g = tr + (1 —t)f* G Bf*,e- The t parameter is chosen at random from the 
uniform distribution on the interval (0,t e ], with t e = e/\\f* — r\\ L 2. 




Figure 8. Schematic representation of our method to con- 
struct deformations of a folding map. First, we construct ran- 
dom maps in the interior of the space of polynomial maps (red 
dots). Then, we use the convexity of PMaps(A n , k) to con- 
struct random convex combinations of interior points and the 
folding map such that the resulting deformations are "close" 
to the folding map (blue dots). 

As we showed in the proof of Theorem 1, PMaps(A n , k) is a bounded 
convex subset of the cone x" =1 /C(A n , A;) C ©f =1 II£. Hence, in order to pa- 
rametrize the interior of the space of polynomial maps we need to first pa- 
rametrize the convex cone of strictly positive polynomials /C_|_(A n , k), whose 
closure is the cone of non-negative polynomials /C(A n ,fc). This is difficult 
to attain in practice, because /C+(A n ,fc) is not finitely generated for values 
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Ph(x) x[J2. 



of k larger than one. However, using Polya's criterion one can construct an 
infinite sequence of finitely generated cones {7Cjv(A n , k)}^ =1 that satisfy 

(B.l) IC N (A n , k) C /Cjv+i(A n , k) C /Qv +2 (A n , k) C ► /C+(A n , fc), 

and converge to /C+(A n , fc) in the limit of large iV. Therefore, in the numer- 
ical applications that concern us here, we use the cone 1Cn(A u , k) with large 
N as an approximation of /C-|_(A n , /c). This can be done more precisely as 
follows. First, we consider the standard parametrization of the vector space 
of polynomials in which every polynomial P(x) is expressed as a linear 
combination of monomials 

(B.2) P{x) = ^ ] Caiet2---a n x \ 1 ' ' ' x n" ■ 

\a\<k 

Then we construct Pjj(x), the homogenization of P(x), which is obtained 
by multiplying every term of degree \a\ < k in Eq. (|B.2j) by the factor 

(Y^1=i ■ Finally, we multiply Ph{%) by (YHH x i) N an d expand the 

resulting polynomial 

fn+l * N 

Xj 

\i=l / 

in the basis of homogenous monomials of degree N + k in n + 1 variables: 

1 X 2 X n x n + l } m=N+k - 

The reader can notice that each term in this expansion corresponds to a 
linear combination of the c Q -coefncients that appear in Eq. (|B.2p . We 
then apply Polya's theorem, which requires each coefficient of the expansion 
to be positive |HLP52j . As the dimension of the space of homogeneous 
polynomials of degree N + k in n + 1 variables is (N + k + n)\/((N + k)\n\), 
imposing that each coefficient in the expansion has to be positive gives rise 
to an equal number of inequalities in n^. We define /Qv(A n ,k) to be the 
semi-algebraic subset that satisfies these inequalities. Furthermore, in order 
to parametrize this convex cone one can determine its generators by means 
of a vertex enumeration algorithm (e.g. [AviOOj ) such that any point in the 
cone can be expressed as a linear combination of the generators with positive 
coefficients. 

In order to learn how one can implement this algorithm, it is useful to 
consider a particular example. Here, we consider deformations of the poly- 
nomial two-fold in PMaps(A 2 , 2). In this case, the space of quadrics has 
dimension six. To determine the plot shown in Fig. 6, we used the value 
N = 8 to construct the cone /Cs(A 2 , 2). The application of Polya's criterion 
gave rise to 66 inequalities in 1R 6 ~ n 2 ,, which define a cone with 900 gener- 
ators (we used the reverse search algorithm Irs to determine the generators 
of the cone [AviOOj ). We scaled each generator by a positive real number 
Aj, such that the maximum value of the corresponding polynomial on the 
simplex was one. Given this basis of scaled generators, we used the Dirichlet 
distribution with 901 concentration parameters 

a\ = «2 = • • • = agoi = 1/1000, 
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to generate random convex combinations of the scaled generators. These 
combinations correspond to rays in /Cg(A 2 ,2). Furthermore, by choosing a 
random number sampled from the uniform distribution on (0, 1) we generate 
random quadrics in the segment of the ray defined by the origin of the cone 
and the random convex combination of the scaled generators. We applied 
this algorithm twice to generate pairs of random quadrics {Rl(x), R^ix)}. 
Then, we computed the maximum of the function R[(x) + R^(x) on A 2 , 
that we denote by S max . It follows that sampling a random number t* 
from the uniform distribution on [0,1/S max ], defines an interior map r € 
PMaps(A 2 , 2) with defining polynomials 

Rx(x) = t*R\(x), R 2 (x) = t*R* 2 {x). 

Finally, the construction of a deformation g of the two-folding map /* re- 
quires an additional convex combination of /* with the map r (see Fig. 8.). 
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